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Abstract. Possible origins of tiie magnetic fields of neutron stars include inheritance from the main sequence progenitor 
and dynamo action at some stage of evolution of progenitor. Inheritance is not sufficient to explain the fields of magnetars. 
Energetic considerations point to differential rotation in the final stages of core collapse process as the most likely source 
of field generation, at least for magnetars. A runaway phase of exponential growth is needed to achieve sufficient field 
amplification during relevant phase of core collapse; it can probably be provided by a some form of magnetorotational 
instability. Once formed in core collapse, the field is in danger of decaying again by magnetic instabilities. The evolution of a 
magnetic field in a newly formed neutron star is discussed, with emphasis on the existence of stable equilibrium configurations 
as end products of this evolution, and the role of magnetic heUcity in their existence. 
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SOURCES OF MAGNETIC FIELD 

At all stages of stellar evolution some percentage of stars 
appears to be strongly magnetic: the magnetic A- and B - 
stars, magnetic white dwarfs, and neutron stars. Do these 
fields have a common origin or are different mechanisms 
at work? What kinds of mechanisms do we know of, and 
which would be the most plausible one for neutron stars? 
Three possible origins have been discussed for pulsar 
magnetic fields: 

• the fossil field hypothesis, also called 'flux conser- 
vation', 

• fields generated by a dynamo process at some stage 
in the evolution of the progenitor star, 

• the thermomagnetic effect in the neutron star crust. 

The last one will not be discussed here, since the effect 
is probably not enough to produce pulsar and magnetar 
fields of the required strength. See, however, Ho et al. 
(2004). 



Fossil fields 

The simplest and most popular hypothesis, flux con- 
servation, is that neutron star magnetic fields are sim- 
ple remnants of the field of their main sequence (MS) 
progenitors. How does this work out quantitatively? The 
strongest fields known in MS stars are around lO^G. For 
a MS progenitor of radius 410"cm, (lOMp)) the star 
would contain a magnetic flux of ~ 5 10^^ G cm^. A neu- 
tron star with the same flux would have surface field 
strength of 5 lO'^G, sufficient for a magnetar. But this 
is an unrealistically optimistic estimate. The neutron star 



contains only some 15% of the progenitor's mass, and 
the inner 1 .4M0 of the progenitor occupy only some 2% 
of the star's crossection. If the field is roughly uniform, 
the flux contained in the core would correspond to lO^^G 
on a neutron star, insufficient for a magnetar. 

The second problem is statistics: only a very small 
fraction of the progenitors has magnetic fields as large as 
lO'^G. Magnetars, on the other hand, are bom frequently. 
As shown in P. Wood's talk (this volume), their birth rate 
is probably comparable to that of normal neutron stars. 

The fossil hypothesis is therefore not very attractive, 
at least not for magnetars. Magnetar fields might have a 
different origin than normal pulsars, of course, but that 
raises the obvious question why a process operating in 
magnetar progenitors could not work, at lower strength, 
in those of pulsars as well. 

In any case, the 'flux conservation' of a fossil field 
would have interesting consequences for the rotation of 
pulsars, since this scenario implies that a strong magnetic 
field is present in the star during its entire evolution. 
This includes the giant/supergiant phases in which the 
envelope rotates very slowly. Since this envelope also 
contains almost all the moment of inertia of the star, any 
couphng between core and envelope will quickly spin the 
core down to the same low rotation rate. The internal 
magnetic field needed to explain the magnetic flux of 
pulsars is more than enough to provide such couphng. 
Except in the very latest stages of evolution, when the 
Alfven travel time through the star becomes longer than 
the evolution time scale, the core will corotate with the 
envelope. The angular momentum left in the core is then 
much too low to explain to rotation rates of pulsars. 

This would disqualify the fossil field hypothesis unless 
the rotation of pulsars is due to something else, for 
example 'birth kicks' , as proposed by Spruit and Phiimey 



(1998). Either the angular momentum, or the magnetic 
field of pulsars can be fossil, but not both. 

Birth kicks would generically lead to correlation be- 
tween spin axis and direction of proper motion (inde- 
pendent of whether the magnetic field is a fossil or 
not). Interestingly, such correlations are now being found 
(Rankin 2006). As in the case of the Crab and Vela pul- 
sars, there appears to be a general preference for align- 
ment of proper motion and spin axis on the plane of the 
sky. In the analysis of Spruit and Phinney, this means that 
the kicks are of relatively long duration, long compared 
with the spin period of the (proto-)neutron star at the time 
of the kicks (Romani 2005, Wang et al. 2006). 



Dynamos 

A magnetic field might be generated in the core of 
a SN progenitor at some earlier stage of its evolution, 
and subsequently amplified during core collapse (for ex- 
ample simply by 'flux conservation'). Possible processes 
are: 



Convection 

Large scale magnetic fields generated by dynamo pro- 
cesses are traditionally associated with convective zones. 
Field strengths of such dynamos are often estimated by 
assuming equipartition of magnetic energy density with 
the kinetic energy density of convection, = A%pv^ 
where v is the convective velocity. The energy flux in 
convection is F « pv^, and if this carries the star's lumi- 
nosity L, the equipartition field strength is found as 

^Mq' 38 8 ' 

where the mean density 3M/AKr' has been assumed for 
p. Compressing such a field from r = 10^ to lO*" cm 
would yield B = lO'^G. This would be sufficient for a 
normal pulsar. 

It is not at all clear that this is a reasonable estimate, 
however In the Sun, the large scale (dipole) field is some 
300 times smaller than equipartition with convection 
would suggest. If this carries over to a stellar interior, 
convection-generated fields in the progenitor would be 
insufficient for producing pulsars fields. 

It is also quite unclear if a field generated by a con- 
vective zone, even if it were strong enough, would leave 
a permanent dipole moment in the core when it becomes 
stably stratified again. Does the magnetic field generated 
in a convective region simply retreat, together with the 
convection, or can it leave a net magnetic field behind? 
Since expansion and retreat of a convective zone is slow 



compared with the characteristic dynamo reversal time, 
it is likely to leave behind, if anything a series of mag- 
netic zones of different orientation instead of a net large 
scale field (much like the famous magnetic stripes along 
midocean ridges in the earth's crust). 



Field generation in stable zones 

Convection is not an essential ingredient in magnetic 
field generation in stars. The energy of differential ro- 
tation is sufficient. 'Closing of the dynamo loop' can 
be achieved by instabilities in the magnetic field itself, 
in much the same way as 'magnetorotational' fields are 
generated in accretion disks (Hawley et al. 1995). 

Contrary to conventional wisdom, field generation can 
thus take place also in stably stratified zones of stars 
(Spruit 2002, Braithwaite and Spruit 2006). The stable 
stratification, however, strongly limits the radial length 
scale on which the process operates (just as it does with 
hydrodynamic processes, cf. Zahn 1983, 1992). For this 
reason it is unlikely to produce a large scale field. 

As in the case of convection it is unlikely that such 
a process would produce a net dipole moment of much 
significance for the neutron star descendant. It would, 
however, be of critical importance for the evolution of the 
angular momentum distribution in the progenitor, since 
even weak magnetic fields with small radial length scales 
can exert torques that are significant on evolutionary time 
scales (Heger et al. 2005). 



Neutrino-driven convection 

The convective velocities in the above estimates are 
driven by the (radiative) luminosity of the star. Much 
more powerful convection takes place during core col- 
lapse, driven by the much higher neutrino flux. The lu- 
minosity is then of order 10^^ erg/s, the size of the core 
some 3 10^, which yields (from eq.[T} I.e. m 

the magnetar range (Thompson & Duncan 1993). How- 
ever, if the actual dipole field generated is as small, com- 
pared with this equipartition number, as it is in the Sun, it 
would be 300 times smaller. Even after contraction from 
3 10^ to lO^cm, this would produce a field of only 10^^ 
G, marginal for a magnetar, though sufficient for ordi- 
nary pulsars. 

Magnetic fields during collapse: energy 
estimates 

The arguments above suggest that field generation in 
the progenitor followed by compression during core col- 



lapse is not the most promising scenario for the produc- 
tion of pulsar fields. The field generation processes in the 
progenitor are important for providing a 'seed field', but 
it still has to be amplified by a large factor 

Suppose a useful initial magnetic field has somehow 
been generated in the pre-SN core. If the collapse pro- 
ceeds with flux conservation, B ^ 1 /R^, the magnetic en- 
ergy Eb ~ B^E? increases as 1 /R, i.e. as a constant frac- 
tion of the gravitational binding energy Eq of the core. 
Eg ^ GM^ /R. Simple flux conservation therefore can- 
not increases a dynamically insignificant initial field into 
a dynamically significant field. A field of of 10^^ G in 
the newly formed neutron star (R = 10^) would require a 
field of 10"G in a pre-coUapse core of R= 10 cm. As I 
have argued above, neither a fossil field nor an internally 
generated (pre-collapse) field is likely to be this high. 

It is possible to do better, at least in simple energetic 
terms, by exploiting differential rotation of the star Un- 
der angular momentum conservation the rotational en- 
ergy varies as 



Ea=\m' = \j^/I-l/R^. 



(2) 



This increases more rapidly than the gravitational energy. 
For the generation of a magnetic field, only the energy 
fraction in differential rotation is relevant: 



EAn = i^) En ~ 0.1£n,~ l/R^ 



(3) 



where I have assumed AH w Q.3Q. as a plausible rate of 
differential rotation. 

Assuming there is a process that can convert all en- 
ergy of differential rotation into magnetic energy Eg = 



^^^NS ^ ^An, this would produce a maximum field 
strength inside the star 
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(4) 



where Hns is the spin of the neutron star end product 
and i^K the Kepler frequency at the neutron star surface. 
The observable surface field strength would be lower 
than this internal field, and the conversion of differential 
rotation into magnetic energy less than 100% efficient. 
These efficiency factors are hard to guess, but it suggests 
that field strengths of the order of a magnetar field are 
plausible, provided the collapsing core contains enough 
angular momentum. 

How large can the angular momentum of the collaps- 
ing core be? An upper limit is set by the average energy 
of observed core collapse supemovae, of the order 10^^ 
erg. If the rotational energy of the neutron star formed 
is larger than this, and the the star strongly magnetic, 
spindown by pulsar emission will be so fast that the ro- 
tational energy is dumped already into the supernova it- 
self This sets an upper limit of 10^' erg to the rotational 



energy, corresponding to a rotation period of 4ms for a 
neutron star of 1 .4Mfj^, or Hns / ^^k ~ 0. 1 . The maximum 
field strength from (Wli is then 3 lO'^ G, assuming again 
AQ./D. = 0.3. This shows that magnetar field strengths, 
up to lO'^, are possible, but require a fairly efficient con- 
version of differential rotation into magnetic energy. 

A way of estimating the angular momentum avail- 
able during core collapse is to follow the evolution of 
the internal rotation in the progenitor up to core col- 
lapse. This requires faith in one's quantitative under- 
standing of the angular momentum transport processes 
in the star. Internal torques due to magnetic fields are 
probably more important than purely hydrodynamic pro- 
cesses alone. Heger et al. (2005) have computed the ini- 
tial spins of neutron stars with the prescription for mag- 
netic torques from Spruit (2002). Rotation periods in the 
range 8-lOms were obtained, somewhat independent of 
the initial (main sequence) spin. This is comfortably be- 
low the upper limit given above, but still rapid enough 
for magnetar fields of lO'^G. 



Consequences and extrapolations 

If differential rotation is an essential ingredient for 
the production of a strong magnetic field, as suggested 
above, at least magnetars must be formed spinning 
rapidly. If a magnetar with B = 3 10 '^G is formed in 
this way, and becomes visible through the supernova de- 
bris after a couple of years, say, it would spin at about 
f ^ 0.1s, and would be more luminous than the Crab 
pulsar It is not clear if this is compatible with observa- 
tions of supernovae in nearby galaxies. 

On the other hand, it is clear from pulsars with inde- 
pendent age estimates (independent from the spindown 
time scale) that normal pulsars are formed with a range 
of initial spins, from ^ 15ms (Crab) to 400 ms (Camilo 
et al. 2007, Gotthelf & Halpern 2007, see also Gotthelf, 
this volume). This shows that there is either a significant 
spread in angular momentum of the collapsing core, or 
that an efficient post-collapse spindown process operates, 
such as friction against a fall-back disk (cf. Li & Jiang 
2007, Wang et al. 2007). 

Some of the slow-born pulsars have rather weak fields 
(Camilo et al. 2007, Gotthelf & Halpern 2007). This 
would agree with angular momentum of the pre-SN 
core playing an important role in determining the field 
strength of a neutron star, as suggested above. 

The rapid initial spin required for magnetar fields can 
be extrapolated to higher rotation rates. At sufficiently 
rapid rotation, the spindown of the nascent, highly mag- 
netized neutron star could power a Gamma-ray burst. 
This magnetar-GRB scenario has recently become popu- 
lar, cf. Kommissarov & Barkov (2007), Bucciantini et al. 
(2007), Yu & Dai (2007) and references therein. 



REQUIREMENTS ON THE 
FIELD-PRODUCTION PROCESS 




FIGURE 1. 'Winding up' of a poloidal field by differential 
rotation. The field lines are stretched into a stronger azimuthal 
field inside the star. The surface field remains unchanged. 

To explain the strength of the dipole component of a 
pulsar magnetic field, not just any field generated inside 
in a SN core will do. An illustrative example is the 
'winding-up' of an initial magnetic field (or 'seed field') 
by differential rotation. 

Suppose the field is initially a purely poloidal field 
Bp (field lines in meridional planes, no azimuthal field 
component), as sketched in fig [T] Differential rotation 
stretches the field lines in azimuthal direction, creating 
an azimuthal field component that increases linearly 
with time, ^ BpAClt. The surface field, and its dipole 
component remain unchanged in this process. 



Stratification 

To produce an interesting pulsar dipole, something 
else has to happen: the strong internal field Bp somehow 
has to be brought to the surface. The obvious process 
is magnetic buoyancy: a strong field reduces the gas 
pressure and density in it, so that loops of the azimuthal 
field tend to float upward against the direction of gravity 
(see fig|2|. The stretch of the loop that finds itself outside 
the star forms a vacuum field with a dipole moment. In 
this way an initially weak poloidal field can be amplified 
by differential rotation into a field with a strong dipole. 

In a stable stratification, a magnetic field can float to 
the surface only if it is strong enough to overcome the 
stratification. This is the case, approximately, when the 





FIGURE 2. Magnetic buoyancy of an azimuthal field loop 
causes part of it to rise, provided the field is strong enough to 
overcome the stratification. This creates a dipole component at 
the surface. 



Alfven speed Va exceeds HN, where H the pressure scale 
height and the buoyancy frequency of the stratification. 
Once the neutron star has formed, the stratification of the 
neutron to proton density ratio is able to prevent fields of 
up to 10^^ G from emerging by such a direct buoyancy 
instability. 

A magnetic field of 10^^ G produced internally by dif- 
ferential rotation must therefore have reached the surface 
in an early stage of the formation of the neutron star, 
when the buoyancy frequency was a percent or less of the 
final value. Buoyant rise is aided, however, by the pres- 
ence of a dense neutrino field (see Thompson & Murray 
2001 for a detailed discussion). 



Instability of poloidal fields 




FIGURE 3. Instability of a poloidal field in a star. (From 
Flowers andRuderman, 1977) 

A poloidal field, such as that created by the buoyant 
process of fig. |2j cannot be the dipole field seen in pul- 
sars, since it is highly unstable. All purely poloidal fields 
in a star are unstable (Wright, 1973, Markey and Tayler 
1974, Flowers & Ruderman 1977). A simple example is 
that of Flowers and Ruderman, shown in fig [3] In this 
example a uniform field inside the star connects to a vac- 
uum field outside. This external field is that of a point 
dipole centered on the star (cf. Jackson E&M). Cutting 
the star in two halves without cutting field lines is possi- 
ble along the axis of the field. Rotating one half in place 
by 180° does not change any of the energy components 
in the interior of the star (thermal, gravitational, mag- 
netic). But the external magnetic energy changes, in the 
same way the energy of two parallel bar magnets changes 
when one is rotated by 180°. The instability proceeds on 
the time it takes an Alfven wave to cross the star 



30 BO 90 120 150 150 210 240 270 500 550 3B0 



50 60 90 120 150 1B0 210 240 270 300 350 350 




How would this bar magnet process work in practice, 
in a gravitating fluid sphere Hke a star? This can be an- 
swered in some detail with numerical simulations. Fig.|4] 
shows the result of calculations by Jonathan Braithwaite. 
The star does not flip hemispheres as fig.[3]suggests. In- 
stead, the instability sets in at a smaller azimuthal length 
scale. The evolution of the surface field distribution re- 
sembles the mushroom-like mixing patterns observed in 
Rayleigh-Taylor instability. Both properties are charac- 
teristic of interchange type instabilities. [The radial field 
component at the surface shows the instability nicely, but 
of course the entire interior of the star takes part in it.] 

In these simulations, the field decayed completely, 
within the numerical uncertainties. This can be under- 
stood from the properties of magnetic helicity, discussed 
below. The conclusion is that buoyant rise of a loop of 
strong field is not enough, at least not in the simple form 
of fig [2] At least one more ingredient is needed to end 
up with a magnetic field configuration that is useful for 
explaining magnetar fields. 



The role of magnetic helicity 

The magnetic helicity ,3^ of a field configuration is a 
global quantity, 

^ = jB-AdV, (5) 

where A is a vector potential of B, and the integral is over 
the volume of the magnetic field. [The gauge dependence 
of A is physically significant. It causes helicity to be 
global: B • A has no meaning as a local 'helicity density']. 
In a perfectly conducting fluid with fixed boundary con- 
ditions magnetic helicity is a conserved quantity. In prac- 
tice perfect conduction is not a realistic situation, since 
rapid reconnection can take place even at very high con- 



FIGURE 4. Numerical simulation of Flowers-Ruderman 
instability in a star. Left panels: radial component of the field 
strength on the star's surface, in the initial dipolar state and at 
two times during the nonlinear evolution of the instability. The 
magnetic equator (radial component zero) is shown as a dark 
line (from Braithwaite and Spruit 2006) 




FIGURE 5. Magnetic helicity. A configuration consisting of 
separated loops of untwisted magnetic field (left sketch) has 
zero helicity. A 'linked poloidal toroidal' field made of the 
same loops (right) is helical. 



ductivity, especially in a dynamically evolving field con- 
figuration. Nevertheless, in laboratory experiments helic- 
ity is often observed to be at least approximately con- 
served. 




FIGURE 6. Equilibrium found as end state of the decay 
of a random field with finite helicity in a stably stratified 
star (idealized sketch). The surface field is determined by the 
poloidal field lines feeding through the hole of the torus. On 
their own they would be unstable as in fig[3j they are stabilized 
by the twisted torus surrounding them. (From Braithwaite and 
Nordlund 2006) 

The consequence is the existence of stable equilibrium 
configurations. To the extent that helicity is conserved, a 
non-equilibrium or unstable magnetic field with a finite 
helicity cannot decay completely, since the helicity of 
a vanishing field is zero. An experiment like that of 



fig |4j a field starting with a finite amount of helicity 
therefore decays initially, but eventually settles into a 
stable equilibrium. [In laboratory plasma physics this 
process of settling is called Taylor relaxation]. 

Two factors are thus necessary for a stable magnetic 
equilibrium in a fluid star: a stable stratification (other- 
wise the field would just escape through to the surface 
by its buoyancy) and a finite magnetic helicity. 

Stellar models with such equilibria have been found by 
Braithwaite and Nordlund (2006). They usually consist 
of an approximately axisymmetric torus of twisted field 
lines inside the star, and a bundle of field lines feeding 
through the 'donut hole' and through the surface of the 
star, see fig. |6] The existence of stable configurations of 
this general form in stars has been surmised several times 
in the past (Prendergast 1956, Kamchatnov 1982, Mestel 
1984). It is quite possible that there also exist stable fields 
with more complex configurations. 



Field amplification in core collapse 

To make use of the energy of differential rotation 
(eq. [3]l, a process has to be found that can convert it to 
magnetic energy, and can do so fast enough. Most of the 
rotational energy becomes available only towards the end 
of the collapse, when the stratification begins to develop 
its stable buoyancy and starts to suppress emergence of 
the field to the surface. This leaves a window of perhaps 
a few seconds for an observable surface magnetic field to 
grow from rotational energy. 

To see what this implies, approximate the collapse as 
a two stage process: an instantaneous collapse to neutron 
star size R = 10^ cm with angular momentum conserva- 
tion and flux freezing, followed by a phase of field am- 
plification by differential rotation. Consider first the sim- 
ple winding-up process of fig[T] With a nominal rotation 
period of 10ms, AQ./Q. ~ 0.3, and a time available of 3s, 
the field can be amplified by a factor 100. To reach a field 
strength of lO^^G, the initial field would have to be 10^^, 
corresponding to 10^ in a pre-collapse core of lO^cm. 
This is marginally achievable with differential rotation- 
generated fields (Heger et al. 2005), but does not leave 
much room for inefficiencies. In addition, as shown in 
fig[T] something more has to happen to convert a wound- 
up internal field into something observable (i.e. a field 
with a dipole moment at the surface). 

This something is likely to be an instability affecting 
the magnetic field. If the instability operates on a mag- 
netic time scale R/v\, it will have noticeable effects only 
when this time scale becomes significantly shorter than 
our time window of a few seconds. If such an instability 
increases the poloidal field strength, the azimuthal field 
strength produced by winding-up in the continuing dif- 



ferential rotation will then increases more rapidly than 
before, and the instability will grow faster The net effect 
is a runaway amplification. 

There are three instabilities that can in principle be in- 
volved in such a scenario. Buoyant instability (fig|2| will 
be most effective when the stratification is convective; a 
Tayler instability (displacements nearly along equipoten- 
tial surfaces, cf. Spruit 1999) is also effective in a stable 
stratification, and finally there is magneto-rotational in- 
stability (MRI) which requires only an outward decrease 
of the rotation rate. 

For the present discussion the main distinction is be- 
tween buoyant and Tayler instability on the one hand, 
and MRI on the other The first two grow on a magnetic 
time scale R / va, hence become effective only when the 
field has already been amplified to some extent by Un- 
ear winding-up. MRI, understood as a form of ('magnet- 
ically enabled') shear instability, grows on a differential 
rotation time scale l/Ail, hence would be expected to set 
in immediately. This gives MRI modes a clear advantage 
when time is of the essence. 



Numerical simulations 

An example of such a process is given by the simu- 
lations of Akiyama et al. (2003), Ardeljan et al. (2005), 
Moiseenko et al. (2006). For an instability with a fixed 
growth rate (independent of the field strength) the time 
t needed to reach a given amplitude B from an initial 
field strength Bq scales as f ~ ln(B/Bo)- This is roughly 
consistent with the evidence given by Moiseenko et al. 
(2006) and suggests that some form of MRI is responsi- 
ble. 

The standard case presented by Moiseenko et al. (a 
ratio of magnetic to gravitational energy of 10-^), cor- 
responds to a wildly optimistic initial field strength (~ 
lO^^G at/; = 10^ or lO'^ at/?= lO^cm). But because of 
the logarithmic dependence on Bq the process still works 
with a more realistic initial field 10^ times smaller. It is 
possible that it will work even better in 3 dimensions, be- 
cause of the additional modes of instability available in 
3D. 



The first three minutes 

Suppose a field of pulsar or magnetar strength has 
been formed during core collapse, and that the stratifica- 
tion of the proto-neutron star has become stable enough 
to prevent the field from escaping by simple buoyant rise. 
After this dynamical phase the field is still out of equi- 
librium or at least far from stable. At the same time the 
star is still fluid: a crust forms only after some 100s, and 



initially is probably more ductile and deformable than 
a cooled-down neutron star crust. The evolution of the 
field over the first few minutes is therefore critical for its 
survival as a pulsar or magnetar field. 

The field will decay by whatever magnetic instabili- 
ties are available to it. Without rotation, these instabili- 
ties will operate on time scales of the order of the Alfven 
crossing time through the star, tp^ « 3s Bj^'. Such insta- 
bilities continue to exist in a rapidly rotating star because 
the Coriolis force is perpendicular to the velocity, hence 
does not contribute an additional energy barrier to be 
overcome by instability (like a stable stratification does). 

The Coriolis force restricts the range of modes that are 
unstable when the rotation rate is larger than the Alfven 
frequency l/tp^. But usually there are some modes that 
remain unaffected by the Coriolis force (a trivial example 
being modes with displacements nearly parallel to the 
rotation axis). Most modes suffer only a reduction of 
their growth rate, characteristically by a factor Q,tp^, while 
some can actually become stable (Pitts and Tayler 1986). 

A simple survival estimate can be made by assuming 
that the growth of the relevant unstable modes is reduced 
by the characteristic factor i2fA (Flowers and Ruderman 
1977). Without the helicity constraint, the field would 
then decay with time t such that the growth time of 
the instability is of the order t: t\Q. « t. After 200s, 
this yields as maximum surviveable field strength in the 
absence of a helicity constraint: 

B(200) wSlO^^GP-i/^^ (6) 

where P is the spin period in s. For an initial spin of 
lOms, fields of the order of observed pulsar strength 
could thus survive until the crust forms and saves them 
from further decay. 

Fields higher than ~ 5 10^^, i.e. magnetar fields, can 
survive only if their helicity is sufficient to reach a stable 
equilibrium. In fact, this holds independent of crust for- 
mation, since lO'^G is about the maximum field strength 
believed to be supportable by a crustal lattice. For the 
survivability of magnetar-strength fields, the neutron star 
can be treated as fluid, with or without crust. [But the 
crust of course affects the actual decay of magnetar 
fields, modulating the energy releases into a series of 
starquakes.J 

The effect of rotation on field decay can again be 
tested by numerical simulation. Braithwaite (2007) finds 
that rotation does indeed slow down the decay, but not 
quite in the simple way assumed above. Starting with 
the uniform intemal field of the Flowers-Ruderman ar- 
gument, oriented at some angle to the rotation axis, there 
is first a stage of decay which proceeds at the Alfven 
crossing time, independent of the rotation rate. The im- 
plication is that this phase is due to the modes that are 
not affected by the CorioUs force. The further decay is 



slowed down by rotation, but within numerical uncer- 
tainty no stable equilibria were found, consistent with the 
vanishing helicity of the initial state. 

SUMMARY 

Three popular scenarios for the origin of magnetic fields 
in pulsars and magnetars are i) the fossil field hypothe- 
sis, ii) fields generated internally in the progenitor, and 
iii) field generation during core collapse. Fossil fields are 
ruled out for magnetars (required field strength of the MS 
progenitor too large, number of magnetic MS progeni- 
tors too small). If a fossil field plays a role for normal 
pulsars, magnetic coupling between core and envelope of 
the progenitor will remove the angular momentum from 
the core, producing neutron stars with extremely long ro- 
tation periods. The spin of pulsars must then be due to the 
'birth kicks' that also give rise to their proper motions. 

Internally generated fields, followed by core collapse 
under flux conservation, is unattractive because of the 
rather low field strengths that can plausibly be produced. 
The most plausible source of energy for field genera- 
tion is the differential rotation during core collapse; it 
increases more rapidly than even the gravitational en- 
ergy. This predicts that the most magnetic stars are also 
formed rotating intrinsically rapidly (but suffering rapid 
spindown by a pulsar wind already during the supernova, 
cf. Metzger et al. 2007). 

The most critical factor for the production of a strong 
field is the finite time available in the final collapse phase, 
on the order of seconds. At earlier times in the collapse 
the energy in differential rotation is not yet high enough, 
at later times (when the neutrinos have left) there may 
still be enough energy in differential rotation, but the 
stratification has then become so stable that even fields 
of magnetar strength cannot reach the surface any more. 

Combined with the lowish initial field strengths in 
the pre-coUapse core, this implies the existence of a 
very efficient field amplification process during collapse. 
Linear winding-up in differential rotation is too slow 
except possibly for normal pulsars. Runaway growth 
of a magnetic field is possible by the combination of 
differential rotation with magnetically driven or shear- 
driven instabilities. Magnetorotational instabihty works 
fastest because it is essentially a shear-driven instability 
operating on the short time scale of differential rotation. 
The simulations of Akiyama et al. (2003) and Ardeljan 
et al. (2005) illustrate the feasibiUty of such a process. 

Most simulations of magnetic core collapse (e.g. Shi- 
bata et al. 2006, Obergaulinger et al. 2006, Sawai et al. 
2007, Burrows et al. 2007, see also the more analyti- 
cal considerations in Uzdensky and MacFadyen 2007) 
start with rather optimistic assumptions on the strength 
and configuration of the pre-coUapse field (an aUgned 



dipole of lO'^G at R — 10^ cm for example). Such ini- 
tial states produce interesting results in themselves. But 
to explain magnetar fields, or magnetically powered su- 
pernovae, let alone the magnetic field strengths needed 
to produce magnetically powered GRB jets, much more 
has to happen than the flux freezing collapse of a dipole 
field. It is likely that pre-collapse fields with more realis- 
tic configurations and strengths, and processes like those 
seen in Moiseenko et al. (2006), will produce magnetic 
fields that behave rather differently. In particular it is still 
an open question how the ordered, rotation-aligned mag- 
netic field configurations are formed that work so well in 
producing jets (or if they form at all). 

Since the field of the neutron star is formed in a highly 
dynamic event, it is likely to be far from equilibrium, 
or at least far from stable. For the same reason it is 
also likely to be dominated by higher multipoles rather 
than the dipole component that determines pulsar phe- 
nomenology. For field strengths of lO'^ G and less mag- 
netic instabilities are slow enough for the field to survive 
until it can be anchored by the newly formed crust (rapid 
rotation helps here, since the Coriolis force reduces the 
growth rate of instabilities). 

Fields of magnetar strength decay rapidly, and are 
too strong to be anchored by the crust. Their apparent 
survival implies that they have found a stable equilib- 
rium configuration within seconds after their formation. 
The formation of such stable equilibria from complex, 
strongly unstable initial fields can be demonstrated with 
numerical simulations. The key ingredients for the exis- 
tence of such equilibria are i) the stable stratification of a 
neutron star, and ii) a nonvanishing magnetic helicity. 
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